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ABSTRACT 

Weak gravitational lensing observations probe the spectrum and evolution of density fluctuations and 
the cosmological parameters which govern them but are currently limited to small fields and subject to 
selection biases. We show how the expected signal from large-scale structure arises from the contributions 
from and correlations between individual halos. We determine the convergence power spectrum as a 
function of the maximum halo mass and so provide the means to interpret results from surveys that lack 
high mass halos either through selection criteria or small fields. Since shot noise from rare massive halos 
is mainly responsible for the sample variance below 10', our method should aid our ability to extract 
cosmological information from small fields. 

Subject headings: cosmology: theory — large scale structure of universe — gravitational lensing 



1. INTRODUCTION 
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Weak gravitational lensing of faint galaxies probes the 
distribution of matter along the line of sight. Lensing 
by large-scale structure (LSS) induces correlati on in the 
galaxy ellipti citics at the p ercent level (e.g., Miralda- 



Escudej 1991; |Blandford ct~al| 1991; |Kaiser| 1992). Though 



challenging to measure, these correlations provide impor 
tant cosmological information that is complementary to 
that supplied by the cos mic microwave back ground and 
poten tially a s precis e (e.g., Jain fc Seljalj 199 7; Bcrnardcau| 
et all 199 7; |Kaiserj 1998 ; ISchneider et al| 1998; |Hu ' 
Te gmarkl 1999; |Coooray| 199 9; |Van Waerbeke et al| 1999; 
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see Bartelmann & Schneider 2000 for a recent review) . In- 
deed several recent studies have provided the first clear 
evidence for weak le nsing in so-called blank fields (e.g 



Van Waerbeke et al 2000; Bacon et al 2000; Wittman et 



; g 2000). 

Weak lensing surveys are currently limited to small fields 
5^ ■ which may not be representative of the universe as a whole, 
| owing to sample variance. In particular, rare massive ob- 
jects can contribute strongly to the mean power in the 
shear or convergence but not be present in the observed 
fields. The problem is compounded if one chooses blank 
fields subject to the condition that they do not contain 
known clusters of galaxies. Our objective in this Letter 
is to quantify these effects and to understand what frac- 
tion of the total convergence power spectrum should arise 
from lensing by individual massive clusters as a function 
of scale. 

In the context of standard cold dark matter (CDM) 
models for structure formation, the dark matter halos that 
are responsible for lensing have properties that have been 
intensely studied by numerical simulations. In particular, 
analytic scalings and fits now exist for the abundance, pro- 
file, and correlations of halos of a given mass. We show 
how the convergence power spectrum predicted in these 
models can be constructed from these halo properties. 
The criti cal ingredients are: the Press-Schechter formal- 
ism (PS; Press & Schcchtci 1974) for the mass function; 
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the NFW profile of Navarro et al (1996), and th e halo bias 



model of Mo fc White (1996). Following 



Scljak (2000), we 



modify halo profile parameters, specifically concentration, 
so that halos account for the full non-linear dark matter 
power spectrum and generalize his treatment to be appli- 
cable through all redshifts relevant to current galaxy el- 
lipticity measurements of LSS lensing. This calculational 
method allows us to determine the contributions to the 
convergence power spectrum of halos of a given mass. 

Throughout this paper, we will take ACDM as our fidu- 
cial cosmology with parameters ft c — 0.30 for the CDM 
density, flf, — 0.05 for the baryon density, fl\ = 0.65 for 
the cosmological constant, h — 0.65 for the dimensionless 
Hubble constant and a scale invariant spectrum of primor- 
dial fluctuatio ns, normalized to galaxy cluster abundances 
(as = .9 see Viana fc Liddle 1999) and consistent with 
COBE ( Bunn fc White| 1997^ For the linear power spec- 



trum, w e take the fitting formula for the transfer function 
given m [Eisenstein fc Hrj (1999). 



2. LENSING BY HALOS 

2.1. Halo Profile 

We model dark matter halos as NFW profiles with a 
density distribution 



p(r) 



(r/r s )(l+r/r s ) 2 



(1) 



The density profile can be integrated and related to the 
total dark matter mass of the halo within r v 



M = Airp s rl 



log(l 



1 



(2) 



where the concentration, c, is defined as r v /r s . Choos- 
ing r v as the virial radius of the halo, spherical collapse 
tells us that M — 4nrlA(z)pb /3, wh ere A(z) is the over- 
density of collapse (see e.g. [Henry 2000) and pb is the 
background matter density today. We use comoving co- 
ordinates throughout. By equating these two expressions, 
one can eliminate p s and describe the halo by its mass Al 
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Fig. 1. — Halo model power spectra, (a) Dark matter power spectrum at redshift of 0.5. (b) Convergence power spectrum with z s = 1. The 
sum of the Poisson (dotted line) and correlation (dashed) contribution (solid line) compares well with that predicted by the non-linear power 
spectrum based on the PD fitting function (long-dashed line). In (a), the linear matter density power spectrum is shown with a dot-dashed 
line. 



and concentration c. Finally, we can determine a relation 
between M and c such that halo distribution produces the 
same power as the no n-linear dark matter power spectrum, 
as outlined in 



Seljak (2000). 



2.2. Convergence Power Spectrum 

For lcnsing convergence, we need the projected surface 
mass density, which is the line-of-sight integral of the pro- 
file 

p(r±,r\{)dr\\, (3) 



where rn is the line-of-sight distance and r± is the perpen- 
dicular distance. As in equation (0), the cut off here at 
the virial radius reflects the fact that we only account for 
mass contributions out to r v (sec Bartclmann 1 996 for an 
analytical description when r v — > oo). The convergence 
on the sky k(6) is related to surface mass density through 
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(l + zi)-E(di6), 



(4) 



where the extra factor of (1 + zi) from the familiar expres- 
sion comes from the use of comoving coordinates to define 
densities and distances, e.g. di, d s and di s are the comov- 
ing angular diameter distances from the observer to lens, 
observer to source, and the lens to source, respectively. 
The total convergence power spectrum due to halos, 



C£ ot , can be split into two parts: a Poisson term, C£, 
and a term involving correlations be tween individu al ha- 
los, C%. This split was introduced by Cole fe Kaisei| (f 988) 
to e xamine the power spec trum of the Sunyaev-Zei'dovich 
(SZ; Bunyaev fe Zcrdovich| f 980 ) effect due to galaxy clus- 
ters (see Komatsu fc Kitayama 1999 and references therein 
for more recent applications). 

The Poisson term due to individual halo contributions 
can be written as, 
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where z s is the redshift of background sources, d 2 V/dzdVL 
is the comoving differential volume, and 



ki = 2tt / 9d9K(9)J 
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(6) 



is the 2D Fourier transform of the halo profile. The halo 
mass distribution as a function of redshift [dn(M, z)/dM] 
is determined through the PS formalism. 

Here, we have assumed that all sources are at a sin- 
gle redshift; for a distribution of sources one integrates 
over the normalized background source redshift distribu- 
tion. The minimum, M m ; n , and maximum, M max , masses 
can be varied to study the effects of rare and excluded high 
mass halos. 

The clustering term arises from correlations between ha- 
los of different masses. By assuming that the linear matter 
density power spectrum, P(k,z), is related to the power 
spectrum of halos over the whole mass range via a redshift- 
dependent linear bias term, b(M, z), we can write the cor- 
relation term as 
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Her e we have utilized the Limber approximation (Lim- 



ber 1954) by setting k = l/d h Mo fc White (1996) 



find that the halo bias can be described by b(M, z) 
1 + [v 2 (M, z) - 1]/6 C , where v(M,z) = 5 c /a(M,z) is 
the peak-height threshold, a(M, z) is the rms fluctuation 
within a top-hat filter at the virial radius corresponding to 
mass M, an d 5 C is the threshold overdensity of spherical 



collapse (see Henry 2000 for useful fitting functions). 



3. RESULTS 



Following the approach given in |Seljak| (2000), we first 
test the halo prescription against the full non-linear den- 
sity power spe ctrum found in simulations and fit by Pea-, 
cock fc Dodds| (PD, 1996). In Fig. |l|a, as an example, we 
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Fig. 2. — Lensing convergence as a function of maximum halo mass. In the increasing order, for both Poisson (dotted lines) and total 
contributions (solid lines), the maximum mass is 10 11 , 10 12 , 10 13 , 10 14 and 10 15 Mq. (a) sources at z s = 0.4; (b) z s = 1; (c) z s = 1 with 
mass cut off applied only out to z = 0.3. 



show the comparison at z — 0.5. A good match between calculated follo wing Hu & Tegmark (1999) using the Pea- 



the concentration relation of 3eljak (2000) as 



the tv r o power spectra was achieved by slightly modifying cock fc Dodds (1996) power spectrum for the underlying 



c(M, z) = a(z) 



M 



M*(z) 



(8) 



Here, M+(z) is the non- linear mass scale at which v (M, z) = 
1, while a(z) and b(z) can be considered as adjustable pa- 
rameters. The dark matter power spectrum is well repro- 
duced, to within 20% for 0.0001 < k < 500 Mpc -1 , out to 
a redshift of 1 with the parameters a(z) = 10.3(1 + z) -0 - 3 , 
and b(z) = 0.24(1 + z) -0,3 , which agree with the values 
given by Seljak (2000) for the NFW profile at z = 0. The 
two power spectra differ increasingly with scale at k > 500 
Mpc -1 , but the Peacock and Dodds (1996) power spec- 
trum is not reliable there due to the resolution limit of 
the simulations from which the non-linear power spectrum 
was derived. Note that the above c(M, z) relation is only 
valid for the cosmology used here and for the NFW profile; 
moreover, it should not necessarily give the true mean den- 
sity profile of halos, since other effects not considered in 
our halo prescription, such as halo substructure, would af- 
fect the relation between the dark matter power spectrum 
and the spatial distribution and mean density profiles of 
halos. A detailed study of c(M, z) as generally applied to 
all cosmologies, profile shapes and power spectra is cur- 
rently in progress (Seljak, private communication). 

In general, the behavior of dark matter power spec- 
trum due to halos can be understood in the following way. 
The linear portion of the dark matter power spectrum, 
k < 0.1 Mpc -1 , results from the correlation between indi- 
vidual dark matter hal os and reflects the bias prescription. 
The fitting formulae of [Mo fc Wriite| (1996) adequately de- 
scribes this regime for all redshifts. The mid portion of the 
power spectrum, around k ~ 0.1 — 1 Mpc -1 corresponds 
to the non- linear scale M ~ M*(z), where the Poisson 
and correlated term contribute comparably. At higher A;'s, 
the power aris es main ly from the contributions of individ- 
ual halos (see aeljak 2000 for a discussion of the detailed 
properties of the density and galaxy power spectra due to 
halos) . 

In Fig. ^d, we show the same comparison for the con- 
vergence power spectrum. The LSS power spectrum was 



mass distribution and using the same Limber approxima- 
tion as the correlation calculation presented here. The 
lensing power spectrum due to halos has the same behav- 
ior as the dark matter power spectrum. At large angles 
(I <J 100), the correlations between halos dominate. The 
transition from linear to non-linear is at I ~ 500 where 
halos of mass similar to M*(z) contribute. The Poisson 
contributions start dominating at / > 1000. 

In order to establish the extent to which massive halos 
contribute, we varied the maximum mass of halos, M max , 
in the convergence calculation. The results are shown in 
Fig. g. We use a background source redshift of 1 and 0.4, 
corresponding to deep lensing surveys and to a shallower 
survey such as the ongoing Sloan Digital Sky Survey^]. In 
Fig. Ha, b we exclude masses above a certain threshold at 
all redshifts and in c only for halos below redshift z — 
0.3, reflecting the fact that current observations of galaxy 
clusters are likely to be complete only out to such a low 
redshift. Assuming the latter, we find that a significant 
contribution comes from massive clusters at low redshifts 
(see Figs. |^b & c). Ignoring such masses, say above ~ 10 14 
Mq can lead to a convergence power spectrum which is a 
factor of ~ 2 lower than the total. Note that such a high 
mass cut off affects the Poisson contribution of halos more 
than the correlated contributions and can bias the shape 
not just the amplitude of the power spectrum. 

In Fig. |^a, we show the dependence of C* ot , for several 
/ values. If halos < 10 15 Mq are well represented in a 
survey, then the power spectrum will track the LSS con- 
vergence power spectrum for all I values of interest. The 
surface number density of halos determines how large a 
survey should be to possess a fair sample of halos of a 
given mass. We show this in Fig. ||b as predicted by PS 
formalism for our fiducial cosmological model for halos out 
to (z = 0.3 and z — 1.0). Since the surface number density 
of > 10 15 Mq halos out to a redshift of 0.3 and 1.0 is ~ 
0.03 and 0.08 degree -2 respectively, a survey of order ~ 
30 degree 2 should be sufficient to contain a fair sample of 
the universe for recovery of the full LSS convergence power 

2 http:/ /www. sdss.org 
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spectrum. 

One caveat is that mass cuts may affect the higher mo- 
ments of the convergence differently so that a fair sample 
for a quantity such as skewness will require a different sur- 
vey strategy. From numerical simulations ( White fc Hrj 

1999), we know that S3 = ^k 3 ) / (k 2 ) 2 shows substantial 
sample variance, implying that it may be dominated by 
rare massive halos. When calcula ted with hyper-extended 
perturbation theory (HEPT, see |Huj 1999) but using the 
halo power spectrum instead of non-linear density power 
spectrum, the skewness decreased by a factor of ~ 2.5 to 
3 with a mass cut off at ~ 10 13 Mq. Since it is unclear 
to what extent HEPT ansatz remains valid for the halo 
description, these results should be taken as provisional 
and will be the subject of future study. 

While upcoming wide-field weak lcnsing surveys, such 
as the ME GACAM expe riment at Canada-France-Hawaii 
Telescope ( Boulade et al 1998), and the proposed wide field 
survey by Tyson et al. (2000, private communication) will 
cover areas up to ~ 30 degree 2 or more, the s urveys that 
have been so far published, e.g., Wittman et al (2000), only 
cover at most 4 degree 2 in areas without known clusters. 
The observed convergence in these fields should be biased 
low compared with the mean and vary widely from field to 
field due to sample variance from the Poisson contribution 
of the largest mass halos in the fields, which are mainly 
responsible for the sample variance below 10' (see |Whitej 
& Hu|l999). 

Our results can also be used proactively. If properties 
of the mass distribution such as the maximum mass halo 
in the observed lensing fields are known, say through prior 
optical, X-r ay, SZ or even inter nally in the lensing obser- 
vations (see Kruse & Schneider 1999), one can make a fair 
comparison of the observations to theoretical model pre- 
dictions with a mass cut off in our formalism. Even for 
larger surveys, the identification and extraction of large 
halo contributions can be beneficial: most of the sample 
variance in the fields will be due to rare massive halos. 
A reduction in the sample variance increases the preci- 
sion with which the power spectrum can be measured and 
hence the cosmological parameters upon which it depends. 

We acknowledge useful discussions with Uros Seljak. 
WH is supported by the Keck Foundation and NSF- 
9513835. 
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